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We investigate the effects of including a quasi-topological cubic curvature term to the Gauss- 
Bonnet action to five dimensional Lifshitz gravity. We find that a new set of Lifshitz black hole 
solutions exist that are analogous to those obtained in third-order Lovelock gravity in higher di- 
mensions. No additional matter fields are required to obtain solutions with asymptotic Lifshitz 
behaviour, though we also investigate solutions with matter. Furthermore, we examine black hole 
solutions and their thermodynamics in this situation and find that a negative quasi-topological 
term, just like a positive Gauss-Bonnet term, prevents instabilities in what are ordinarily unstable 
Einsteinian black holes. 



I. INTRODUCTION 



The concept of holography has proven to be enormously fruitful for demonstrating interesting new connections 
between disparate areas of physics. The basic idea is that gravitational dynamics in a given dimensionality can 
be mapped onto some other (non-gravitational) field theory of a lower dimensionality. Holography has been most 
thoroughly explored in the context of the AdS/CFT correspondence conjecture, in which a large volume of calculational 
evidence indicates that a (relativistic) conformal field theory (CFT) can be mapped to gravitational dynamics in an 
asymptotically Anti de Sitter (AdS) spacetime of one larger dimension 

Over the past few years it has become clear that holographic concepts cover a much broader conceptual territory. 
For example holographic renormalization has been shown to be a useful tool for understanding conserved quantities 
and gravitational thermodynamics in both asymptotically de Sitter Q and asymptotically flat spacetimes |3| . Much 
more recently holography has been extended to describe a duality between a broad range of strongly coupled field 
theories and gravity in the context of QCD quark-gluon plasmas [J], atomic physics, and condensed matter physics 

Gravity- gauge duality is evidently a robust concept, and its full implications for physics (for example in elucidating 
the strong coupling behaviour of the non-gravitational theories noted above) remain to be understood. One line of 
investigation has been concerned with Lifshitz field theories, which have an anisotropic scaling of the form 

t -> X z t, r -> A~V, x -> Xx (1) 

exhibited by fixed points governing the behaviour of various condensed matter systems. Such scaling offers some 
promise for further extending holographic duality between condensed matter physics and gravity. While for z = 1 this 
scaling symmetry is the familiar conformal symmetry, for z = 3, theories with this type of scaling are power-counting 
renormalizable, possibly providing a UV completion to the effective gravitational field theory [8|. For Lifshitz gravity, 
the natural asymptotic spacetime metric is 

^ = -^ 2 + L- 2 (^+,w) (2) 

noted earlier in a braneworld context [9(. 

A _D-dimensional anisotropic scale invariant background using an action that couples gravity to a massive gauge 
field (or alternatively to 2-form and dualized (D — l)-form field strengths with a Chern-Simons coupling) can be 
constructed that has solutions with the asymptotic behaviour ([2]) [l(| ■ An early example [ll[ for an extended class of 
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vacuum solutions for a sort of higher-dimensional dilaton gravity with general z was soon followed by the discovery 
of black hole solutions, both exact (for z = 2) [l2[ and numerical (for more general values of z) [i"2T - [l5j ]. 

Since in general one expects quantum-gravitational effects to induce corrections to the Einstein action, it is natural 
to consider modifying the gravitational part of the action with higher-derivative terms due to additional powers 
of the curvature. Such terms must be considered on the gravity side of the duality conjecture in order to study 
CFTs with different values for their central charges. Here, Lovelock gravity theories play a special role in that the 
number of metric derivatives in any field equation is never larger than 2. Furthermore, third-order Lovelock gravity is 
supersymmetric, and therefore one can define superconformal field theories via the AdS/CFT correspondence [16lll7j|. 

The addition to the action of a term cubic in curvature is not new, but asymptotic Lifshitz solutions in Lovelock 
gravity coupled to a massive Abelian gauge field were only recently discovered [18|. For a suitable choice of coupling 
constant one can dispense with this massive gauge field, since the additional Lovelock terms can play the role of the 
desired matter. Some new exact black hole solutions were obtained as well as a broad class of numerical solutions, 
and asymptotic Lifshitz solutions with curvature-squared terms in the action have also been investigated [2p| - [23| . 
Somewhat remarkably, the relationship between the energy density, temperature, and entropy density is unchanged 
from Einsteinian gravity [24| even though the sub-leadin g la rge-r behaviour of Lovelock-Lifshitz black branes differs 
substantively from their Einsteinian Lifshitz counterparts TU, 15 1. The relationship between entropy and temperature 
is also the same as the Einsteinian case, apart from a constant of integration that depends on the Lovelock coefficients. 

The hallmark feature of Lovelock theories is that no field equation has more than two derivatives of any metric 
coefficient. However, a new cubic curvature term was recently introduced that can perhaps be regarded as a general- 
ization of Lovelock gravity in five dimensions [25l [26|. The generality arises from a spherical symmetry requirement: 
the field equations will generally reduce to second-order system of differential equations when the metric is spherically 
symmetric. 

This particular class of correction terms has been coined quasi-topological gravity, since in some ways they behave 
like topological invariants in 6 dimensions, yet for nonspherical geometries, they contribute nontrivially to the action. 
Furthermore, there are no Lagrangians that are cubic in curvature in four dimensions for spherical symmetry that 
lead to second order differential equations. Quasi-topological gravity has been previously studied in the case of planar 
AdS black holes. In this paper we investigate the implications of this new term for asymptotically Lifshitz spacetimes. 

Specifically, we examine the effects of higher-curvature modifications to Einsteinian gravity to asymptotically Lif- 
shitz metrics, both with and without massive background Abelian gauge fields. We find that quasi-topological Lifshitz 
gravity replicates the field equations from third-order Lovelock-Lifshitz gravity [18j . provided the quasi-topological 
parameter fj, is appropriately renormalized. We find that indeed, asymptotic Lifshitz black holes exist in both cases. 
We obtain both exact solutions and numerical ones, the latter obtained via the shooting method. We close with a 
short discussion of the relevant thermodynamics and conserved quantities of our black hole solutions. 



II. QUASI-TOPOLOGICAL GRAVITY 

The quasi-topological additions consist of 3rd-order curvature corrections to Gauss-Bonnet gravity that maintain 
second-order field equations with respect to the metric under conditions of spherical symmetry. We use the action 

/ = / d D xyf~ g (-2A + d + ^ -C 2 - t rC 3 - -F UV F» V - -tti 2 A u aA (3) 

J v y V (£>-3)(D-4) (£>-3)(D-6) 4 2 J w 

where D is the number of dimensions (larger than 4 and different from 6), F^ v = d^A^, /i and A are the correction 
terms' coefficients, C\ = R is the Ricci scalar, £2 = R^v-ysR^ 15 — 4i? M „i? A ' 1 ' + R 2 is the Gauss-Bonnet Lagrangian, 
and £3 is the quasi-topological gravity correction. This quasi-topological gravity correction has the form 

r _ 2-P — 3 n » Cn t »n /i * 1 \ ( ~ 8 p \ vp 

~ 3-D 2 — 15-D + 16 ^ A p a (D — 4)(3-D 2 — 15-D + 16) \ 8 ^ vp 

-{D - 2)R^ p R^\RP T + DR^ vp R^R x " + 2(D - 2)R^R X »R^ - ^£-±R^R^R + Rj^j (4) 

This term is only effective in dimensions greater than 4, and it becomes trivial in 6 dimensions [26| . 

Rather than write down the full tensorial expression for the field equations, as we are interested only in spherically 
symmetric solutions we will insert the Lifshitz metric 

L zz r £ g[r) 
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into the action and then functionally vary it, obtaining (after eliminating redundancies) three equations of motion for 
the two metric functions and the gauge field. Boundary conditions require that f(r) and g(r) asymptotically reach 
unity. The term dil 2 is the metric for a constant curvature hypersurface 



D-2 i-l 

dn 2 = de 2 + k- 1 sin 2 (Vfcfli) ( de 2 2 + sin2 B 3 dBl 

i=3 j=1 



(6) 



where parameter k is either —1, 0, or 1, providing hyperbolic, flat, and spherical geometries, respectively. For k = 
a coordinate transformation will reduce this to the form ^2^~ 2 dOu 2 ■ Symmetry requirements imply that the gauge 
field ansatz is 



(J) 



with all other components vanishing. 

Now that the formalism has been specified, we restrict our considerations to five dimensions and so (unless otherwise 
stated) the following results are only valid D = 5. Rather than carry out a full variational principle, we insert the 
ansatz (JSJ) and ([7]) into the action, obtaining the effective action 



/ d4x J dr ^i7^\~g ({ 3r4 {~T l2 - k + Xk2 + ^ 3 ) }' + S Yf{ 9 ( rh ' + zh ? + m2L " h2 ) 



(8) 



for the spherically symmetric case, where n = yg — 

Functionally varying ([5]) with respect g(r), /(r), and h(r) respectively yields upon simplification 

AL 2 r 6 + (3z + 3) r 6 g - 6z\r 6 g 2 + 6z\r 4 L 2 kg - 3r 4 L 2 k - (9z - 3) fiT 6 g 3 + (18z - 9) 
Hr 4 L 2 kg 2 - (9z - 9) ^iL A k 2 r 2 g - 3 f iL 6 k 3 + g(\nf)' i\ 7 - 3 \r 7 g + 3\r 5 L 2 k - ^r 7 g 2 



9nr*gL 2 k - ^r 3 L 4 k 2 



q 2 r 6 



3r 4 



A 



K + Xk + fiK 



q 2 r 3 



zh 



m 2 L 2 h 2 



2/ 



g[rh + zh) + m 2 L 2 h 2 
h 



2r 2 h" - r [(In/)' - (In g)'] {rti + zh) + 2(z + 4)rti + 6zh = 2m 2 L 
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(9) 
(10) 
(11) 



where a prime (') represents differentiation with respect to the radial coordinate r. 

Before trying to find solutions to the above equations, we present a first integral for the above equations of motion. 
It is a matter of calculation to show that this conserved quantity can be written as 



C = [(1 - 2Xg - 3 M5 2 ) (rf + 2 (z - 1) /) - q 2 (zh + rti) h] 



„z+D-2 



1/2 



(12) 



with details of this result given in appendix 1X1 For z = 1, f(r) = g(r) and the constant reduces (in the matter-free 
case) to 



Co 



which is proportional to the mass of black hole. 



V- 



(/-A/ 2 -m/ 3 )', 



III. BLACK HOLES 



A. Matter-free Solutions 



Setting h(r) = 0, we first consider the existence of solutions of the form 



,,2 z 



ds z 



^-dt 2 

L 2z 



2^2 



(13) 



i=i 
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where k = 0. This is a Lifshitz analogue of flat space for asymptotically flat solutions, whose properties have been 
discussed elsewhere [27]. We shall refer to such solutions as "Lifshitz solutions". 
For the metric (fl~3| the field equations (J9j [101) imply 

A = -|j(2 -A), M = i(l-2A), (14) 
independent of our choice of z. These constraints reduce to those of five dimensional Gauss-Bonnet gravity [18j 

K =-J2 aild X =\ ( 15 ) 

when fi — 0. Note that the same constraints are necessary to ensure the existence of asymptotic Lifshitz solutions if 
k £ 0. 

With the above constraints, the exact Lifshitz solution (that is, f(r) — g(r) = 1) is a solution to the field equations 
for any value of z. Eq. (|10p with the condition (TT4")) reduces to 

C 

2- A-3k + 3Ak 2 + (1-2AW 3 = — (16) 
where C is a constant of integration. For C = 0, n = 1 or 

sM = i + ^r ( 17 ) 

yielding the only solution of eq. (fT6|) that has the desired asymptotic behaviour. The function f(r) is not restricted 
by Eq. ©. This degeneracy of the field equations has been noted previously in 5 dimensional Einstein-Gauss-Bonnet 
gravity with a cosmological constant [lj| and third order Lovelock gravity [l~8|. In the Gauss-Bonnet case, it was 
shown that there exists a degenerate set of solutions where f(r) is left unspecified, while for certain values of the 
Gauss-Bonnet parameter, f(r) = g(r). In our case, this degeneracy is lifted when matter is present, and we obtain a 
family of solutions that become unique for a specific field strength, as we shall see. 

Choosing f(r) = g(r) (as in Lovelock gravity |18l |) yields for k = — 1 an event horizon, and consequently the metric 

r 2z / t2\ T 2 rlr 2 

d s 2 = - r —[l-± T )dt 2 + Ldr T2 +r 2 dn\ (18) 
T 2z \ r 2 I -2M - ^ _1 



(1- 



which is an exact black hole solution. 

For C 5^ 0, one can find K and therefore g(r), but upon inserting this solution in eq. ©, we find that the solution 
for f(r) does not exhibit the desired asymptotic behaviour. We find, with one exception, no other exact solutions to 
the field equations for these symmetries and asymptotic behaviour. 

The exception is z = 1 (AdS), for which an exact solution can be found. The requirements that /(r) = g(r) and 
h(r) — produce exact solutions dependent on A if \l = [30j |. Setting \x ^ 0, we first seek solutions for z = 1 without 
any background gauge field. Restricting f(r) = g{r) and setting h(r) = 0, the field equation (fTTj) disappears, while 
the equations (O and (fTTTj) are not independent and can be analytically solved. The result is 



f(r) - g(r) = — - — + 

r 2 3/i 12/ir 2 

where we define 



(Vr + J2(r) + J(r)) 3 - (Vr + J 2 (r) - J(r) 



(19) 



r = -(l6r 4 (3^ + A 2 )) J J(r) = 16r 6 (^4A 3 + 18/xA - 9fi 2 AL 2 - 18^^ J (20) 

and M is a constant of integration. This solution matches the form of one obtained in 3rd order Lovelock gravity for 
D > 6 [30] , as our field equations are of the same form. With this exact solution, we are able to compare results with 
the numerical algorithm. 



B. Matter solutions 



In the presence of a massive gauge field (h(r) ^ 0) the Lifshitz solution (fT3f is also supported by quasi-topological 
gravity provided 

2 2(z-l)(l-2A-3u) 2 3z 

q= z m =T 2 

A = --^ [(l-2A-3Ai)(2z + z 2 )+9-6A-3/i] A<i(l-3/x) (21) 
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where the last constraint arises because we require q 2 > 0. This in turn implies 

3 M (z + l) 2 -(z 2 + 2z + 9) 3 
2^ <A<- — (1 + M ) 

provided A > 0, as is normative for Gauss-Bonnet gravity in the context of heterotic string theory [32j |. 

We look for black hole solutions using both near-horizon and asymptotic series expansions of the metric and gauge 
functions. The near-horizon series solutions are then used to obtain initial conditions for the numerical solution of the 
field equations. The restrictions on /(r) and g(r) now become more rigid: these functions must not only approach 
unity as r — > oo (to satisfy the asymptotically Lifshitz boundary conditions) but they must also tend towards zero 
as r — > r in order to ensure an event horizon exists. First, we will show that series representations exist near and 
far from the horizon, and then we present a set of solutions obtained by numerically solving the differential equations 

nam). 

1. Series Solutions 

We begin by searching for well-behaved black hole solutions in a near-horizon regime. Our ansatz requires that the 
metric functions go to zero linearly near the horizon r = tq: 

f(r) = fi{(r- r ) + /a(r - r f + f s (r - r ) 3 + ...} , 

9(r) = Si (r - r ) + g 2 (r - r ) 2 + g 3 (r - r ) 3 + (22) 
h(r) = fl' 2 {h + h 1 (r- r„) + h 2 (r - r Q ) 2 + h 3 (r - r ) 3 + ...}, 

and we find that the substitution of this ansatz into our equations of motion results in ho — and a restriction on gi : 



z 
ro 

+2Ar 6 ((z + l) 2 + 2) - r 4 (r 2 {z + l) 2 + 8r 2 + 6L 2 k) } 
[3/i ((1 - z)/i! 2 r 5 - 3L 4 k 2 z) + 2A ((1 - z)^ 5 ^ 2 

+3zL 2 kr 2 ) + ((z - l)/n 2 r 5 + 3zr 4 )] ^ . (23) 

We are left with two free parameters, hi and f±, and values for these are selected to ensure proper asymptotic 
behaviour for large r. All of the other terms in the expansion are solvable in terms of these two parameters. 
Solutions at large r can be obtained by linearizing the system, using the ansatz 

/(r) = l + e/ e (r), 
g(r) = l + eg e (r), 

h{r) = l + eh e (r), (24) 
yielding rather lengthy expressions for the leading terms. We have relegated these to Appendix [Cl 

2. Numerical Solutions 
We will find it easier to obtain numerical solutions by writing 

f (25, 



91 =— {3li (r„ 6 ( 2 - l) 2 + W - 2LV 
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-^—l 12 ^ 6 /;? 3 + \-6L 2 kn + 4Ar 2 l fg 2 



in which case the set of differential equations (l9llTTj) can be written as 

dj _ zh(3- 2g) - z 2 hg - rjg (2z + 3) r 2 h 2 (zh + rj) (z - 1) 
dr gr 2 fgH 

fr = ^{^ + ^- 1 ^^ 
+ 3zr 4 [-3fc^i 2 (4z - 2) + (4z)Ar 2 ] fg 2 
~ 3zr 2 [3(2z - 2)fc 2 /i + 4zfcAL 2 r 2 + (2z + 2)r 4 ] /g 

- zfi {[3(z - f) 2 + 12z]r 6 - 6fc 3 L 6 } / 

- zAr 2 { [2(z - l) 2 + 8z + 4]r 4 } / 
+ zr 4 {[(z - l) 2 + 4z + 8]r 2 + 6fcL 2 } / 

+ (z- l)L a r 6 [(zh + rjfg - 3zh 2 ] j 
d£ _ 

dr 3zr 3 fH 

- 3zr 2 [Ak\L 2 r 2 + 4r 4 ] fg 

- zh{[3(z - I) 2 + 12z]r 6 - 6fc 3 L 6 } / 
-zAr 2 {[2(z- l) 2 + 8z + 4]r 4 }/ 

+ zr 4 {[(z - l) 2 + 4z + 8]r 2 + 6kL 2 } f 

+ (z- l)L r e [(zh + rj) 2 g + 3zh 2 ] } , (26) 

where for simplicity we define L = — 1 + 2A + 3/i and H = r 4 + 2\r 2 (kL 2 — r 2 <?) — 3/i (kL 2 — r 2 gy . 

Equation (|25H26p form a system of four coupled first order ordinary differential equations. With these ODEs, initial 
conditions are chosen from the series solution (evaluated just beyond the horizon), and then the shooting method 
(explained in [T3. fl5l|) is used to obtain solutions. 

We consider values of fi and A that guarantee positivity of the energy flux in the dual conformal field theory [Hj when 
z = 1. For z =/= I the dual theory is not well understood and the analogous allowed ranges of \i and A are not known. 
Furthermore microscopic constraints such as positivity of energy and causality are not necessarily responsible for 
setting the lower bound on the ratio of shear viscosity to entropy density in the plasma, since hydrodynamic transport 
is determined by the infrared properties of the system, which do not necessarily enter into the microcausality analysis 
of the theory [31] . However for the most part we shall employ the same values of /x and A as for the z — 1 case, noting 
departures from these values for illustrative purposes as appropriate. 

The specific case z = 1 eliminates the charge q, and the solution is given by equation (fT9|) with f(r) = g(r). Solving 
the system (f9lfTTj) yields numerical solutions. 

We can check the validity of our numerical approach by comparing this to the exact solution in equation (1191) . For 
example, for ro = 0.9, and k = we see from Figure[T]that the two curves (numerical and analytic) are coincident. To 
be certain, we tested equality of the two approaches for /j, = —0.001, A = 0.04, k = 0, ro = 1.5. Evaluating between 
r = 1.51 and r = 15 at intervals of 0.01, we find that the two solutions differ by no more than 10 -7 . 

For z 7^ 1, we numerically obtain solutions for large, medium, and small values of ro over a broad range of initial 
values of the field strength (hi). The quantity f\ is then fixed by asymptotic conditions. For a given value of hi, we 
find that large black holes are asymptotic to functions that monotonically tend to unity, whereas the metric functions 
for small black holes exhibit a spike in magnitude before settling down. However, due to the extra degree of freedom 
in the gauge field strength, we can obtain a family of solutions (and control the spike) by varying hi , subsequently 
adjusting fi to satisfy the asymptotic conditions. In Figure [H we see the result of varying the initial value of hi from 
2.6 (dashed solution) to 2.8 (solid solution). For these the initial values of fi remained constant at 2.0. Note that the 
initial spike present for hi = 2.6 vanishes for hi = 2.8. 

In Figure[3l we plot the metric and gauge functions for a large black hole. All three functions monotonically increase 
from zero at the horizon to unity for large r. Figure |4] shows a medium black hole (ro = 2.4), where the dashed line 
is Einsteinian gravity and the dotted line is quasi-topological gravity. Here, due to the more favourable scale, we see 
that the solution for h(r) is noticeably different. The scale is still too large to see any effect on the g(r) solution, 
however. 

We can see from Figure [S] that for small black holes, f(r) spikes sharply. The plot shows a comparison between 
Einsteinian gravity (dashed), Gauss-Bonnet gravity (solid), and quasi-topological gravity (dotted) for k — — 1. Small 
black holes for k = and k — 1 exhibit similar behaviour. 
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f(r). 



FIG. 1: Comparison of Analytic versus Numerical Solution for f(r), where ro = 0.9, A = .04, /j, = —.001, and k = 0. The two 
curves are identical to one part in 10 -7 . 



The plot in Figure [5] better shows the effect of larger values of \x and A, elucidating how the quasi-topological term 
really affects solutions. 




1 2 3 4 f 5 fi 7 8 9 



FIG. 2: Comparison of two sets of z — 2 solutions for fei = 2.6 (dashed) and hi = 2.8 (solid) for A = .1 and /j, = .001, where 
f{r),g(r),h(r) are plotted versus r respectively in magenta, green and blue. 



IV. THERMODYNAMICS 



In this section we generalize from 5 to D dimensions to study the thermodynamic behaviour of the solutions we 
obtain. The Iyer/Wald prescription for black hole entropy is [29[ 

S=-2n<f d D - 2 x^Y abcd e ab e cd , where Y abcd = -|^— (27) 

J Otiabcd 

where e a b is the binormal to the horizon and C is the Lagrangian, with Latin indices denoting quantities projected 
onto the horizon surface. For the static black holes considered here, Y = Y abcd i a te c d is constant on the horizon and 
so the entropy is given simply as 

S = -2ttY J d D - 2 xy/§, (28) 

where the integration is done on the (D — 2)-dimensional spacelike hypersurface of the Killing horizon with induced 
metric g a b (whose determinant is g). Although the asymptotic behaviour of our solution is different from Ref. (26j 
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FIG. 3: Large black hole (k = 1), where A = .04 and fi — —.001, with f{r),g(r),h(r) versus r respectively in magenta, green 
and blue for z = 2. 



f&) 



h(rj 



/./ 



FIG. 4: Medium black hole of radius ro = 2.4 (k = 0); here A = .04 and /i = —.001, and f(r),g(r),h(r) versus r for z — 2 
respectively in magenta, green and blue. The dotted line is Einsteinian gravity and the solid line is quasi-topological gravity. 



and f(r) ^ g(r), we obtain the same result: 



A f< ■ 2 ^- 2 ) Afc g- 3 ^-y^) (29) 



AG D V D-4 r 2 Q D-6 

where D is the number of dimensions and A is the surface area of the black hole (since our metric is spherically 
symmetric, the surface area will be proportional to r ~ 2 ). 

The temperature of the black holes is found by ensuring regularity at the horizon after Wick-rotation; we obtain 

r -(^L- < 3o » 

This quantity can be numerically calculated, and plotted against entropy on a logarithmic scale, to study stability of 
the black holes. A negative slope indicates that the black hole will not be in thermal equilibrium and so must decay. 
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FIG. 5: Small black hole ro = 0.9 (k = —1), where A = .04 and n — —.001, with f(r),g(r),h(r) versus r for z — 2 respectively in 
magenta, green and blue for Einsteinian gravity (dashed), Gauss-Bonnet gravity (solid), and quasi-topological gravity (dotted). 

/-\ 
- / \ 
i \ 
! \ 

5-1 \ 




4 6 s 10 r 12 14 16 IS 20 



FIG. 6: Medium z — 2 black hole (ro = 2) with /j, = 2.5 and A = — 10 for k = — 1 where Einsteinian gravity is solid, Gauss- 
Bonnet gravity is dashed, and quasi-topological gravity is dotted. Similar to above, f(r0,g(r),h(r) are magenta, green, and 
blue. 

A. Stability of AdS Black Holes 

Plotting the solution for z = 1 in five dimensions, which can be checked with the analytic case, we obtain Figure 
[7J where we use A = 0.4 and /i = —0.001. The solid line is quasi-topological gravity, while dots correspond to Gauss- 
Bonnet and crosses are Einsteinian. The parameter k varies between —1, 0, 1, coloured green, blue, and magenta, 
respectively. Up to the black hole sizes for which we are able to find valid numerical solutions, we see no evidence of 
unstable black holes for any value of k. In the Einsteinian case we see that small black holes will become unstable for 
k = 1, so it is expected that for sufficiently small values of A and fi, the solution will be one of small, unstable black 
holes. 

To see what effect the sign of the quasi-topological parameter has on black hole stability, we plotted a similar set 
of curves for a positive value of /i = 0.001. This plot is shown in FigurelU We sec that for k = +1 sufficiently small 
black holes are thermodynamically unstable. 

B. Stability of Lifshitz Black Holes 



For z = 2, we also plot log(T) versus log(5). Note that this plot is also specific to the five dimensional case. For 
large black holes, it appears that by varying z we do not change the temperature-entropy relationship, but instead 
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FIG. 7: log(T) versus log(S') for z = 1 , A = .04, and /x = —.001. The solid line is quasi-topological gravity, the dotted Gauss- 
Bonnet, the dashed Einsteinian. The parameter k varies between —1, 0, 1, coloured turquoise, red, and brown, respectively. 




FIG. 8: log(T) versus log(.S') for z = 1, A = .04, and /i = 0.001. The parameter k varies between —1, 0, 1, coloured green, blue, 
and magenta, respectively. 

merely introduce a scaling factor to both entropy and temperature terms. It is also apparent that in both cases, 
positive Gauss-Bonnet and quasi-topological terms will both introduce k = +1 instability in black hole solutions. A 
sufficiently negative quasi-topological term is also seen to partly counteract the positive Gauss-Bonnet term in Figure 

M 

V. CONCLUSIONS 

It is well-known that the third order Lovelock term (cubic in the Riemann tensor) does not appear in the field 
equations in five dimensions as it is a topological invariant. Terms cubic in curvature in general yield higher-order 
differential equations for metric components. Quasi-topological gravity |26| is an exception to this general rule - the 
cubic terms conspire to yield second-order differential equations for spherically symmetric metrics. 

The main result of our paper is to demonstrate that a broad class of solutions - those that are asymptotic to 
Lifshitz gravity - exist in quasi-topological gravity in five dimensions. We obtain a family of solutions dependent on 
two parameters, one giving a measure of the gauge field strength and the other the black hole radius. For a given 
value of the gauge field strength, we found that there exists a unique solution with asymptotic Lifshitz behaviour. 
Varying the gauge field strength, we found that there exists a family of solutions for a given black hole radius. The 
r-dependence of these solutions varies considerably: the metric functions can develop a "spike" by increasing the 
gauge field strength. We also find that in general, the quasi-topological term acts similarly to the Gauss-Bonnet 
term, but in negative sign. When a negative Gauss-Bonnet term decreases the magnitude of the spike, a positive 
quasi-topological term will have the same effect. We see this when a positive quasi-topological parameter is added 
to a Gauss-Bonnet solution that has decreased the magnitude of an Einsteinian spike in g(r): our quasi-topological 
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FIG. 9: log(T) versus log(S) for z = 2, A = .04, and = —0.001. The parameter k varies between — 1, 0, 1, coloured magenta, 
blue, and brown, respectively. 
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FIG. 10: log(T) versus log(5) for z — 2, A = .04, and [i — —.0003. The parameter k varies between —1, 0, 1, coloured magenta, 
black, and brown, respectively. 

parameter further decreases the magnitude of the spike. 

We also investigated the thermal stability of these quasi-topological Lifshitz black holes. We found that a negative 
quasi-topological term, just like a positive Gauss-Bonnet term, will prevent instabilities in what are ordinarily unstable 
Emsteinian black holes. For the asymptotically AdS case {z = 1) we found that for sufficiently negative values of fi 
the instabilities that arise in Einsteinian gravity may be removed, in the same way that a sufficiently positive Gauss- 
Bonnet term removes the small-r black hole instability. The AdS solutions were seen to be unstable for positive 
values of [i. With regard to the stability of asymptotically Lifshitz solutions with z — 2 in Einstein gravity [Hj], 
we found that the quasi-local term with positive coupling constant furnishes an effect similar to a positively coupled 
Gauss-Bonnet term. 

It is clear that there is much to explore in Lifshitz gravity with the addition of higher-order curvature corrections. 
The implications of these corrections for the dual theory remain an interesting subject for future investigation. 

Appendix A: The conserved quantity along the radial coordinate r 

In this appendix, we demonstrate the existence of a constant Co, which conserved along the radial coordinate r, and 
compute its form. Since there is no exact Quasi-topological-Lifshitz solution (except under special circumstances), we 
calculate it at the horizon and at infinity. 
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Reparametrizing the metric with the relations 

F(r) 
G(r) = 



1 r 
-ln/(r) + zln-, 

1 r 
-\ng(r) - In-, 



the metric becomes 



ds 2 



R(r)=hx~, 

H(r) = In h(r) + z In • 



^Hr) dt 2 + e 2G(r) dr 2 + ^r) J_ ^2 



(Al) 



(A2) 



whose form we insert into the action. Following a similar method to reference [241 ] . we integrate by parts and obtain 
a one-dimensional Lagrangian Lm = C-ig + £im as 



C lg = {D-2) -2 



A 



„2G 



£>-2 



2F J R + (D - 2)i?' 2 



AL 2 



[4F'i?' 3 + (£> - 5)i?' 4 ] e 



Ml „-2G 



- gi 4 [6F'ii /5 + (D- 7)R' 6 ] e- 4G J e ^-G+(o-2)fl 

£ lm = l -q 2 (m 2 + H' 2 e- 2G ) e -F+G+(D~2)R+2H ^ 



(A3) 



We are then able to write the equations of motion in the same manner as (24| . obtaining the conserved quantity 
C = 2 (F' — R') (1 - 2AL 2 i?' 2 e- 2G - 3 M L 4 i?' 4 e - 4G ) e F - G +^ D ^ R 

- q 2 H' ' e - F -G+(D-2)R+2H 



[(1 - 2A.g - 3 R g 2 ) (rf + 2 (z - 1) /) - g 2 (zh + rti) h] 



n z+D-2 



1/2 



(A4) 



This derivation was performed using the form of the quasi-topological Lagrangian for D dimensions, given by ((3]), and 
the form of the conserved quantity was checked explicitly for D from 7 through 11 to determine the dimensionally 
independent form given. For any value of z, this conserved quantity arises from the symmetry 



/ F(r) 
R(r) 
G(r) 

\H(r) 



For z = 1, f(r) = g(r) and the constant reduces to 



( F(r)+5 \ 
R(r) j±- 2 
G(r) 
\ H(r)+5 J 



(A5) 



Co = ^(/-A/ 2 - / ,/ 3 )' 
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Appendix B: Near-horizon Series Solution Coefficients 



Here we write down the remaining coefficients of the near-horizon series solution (I22[) up to second order. Defining 
for simplicity L = — 1 + 2A + 3/i, we obtain 



+ ^hi 2 L^j + 12zr 7 ^z (g! + ]^L Q hi^j + ^gi - ^L hi 



h = ^-6.gi r o 8 (z (\gi - ^h! 2 L 

+18zr 6 Qj L oz 2 + ^L z + ^L + fiL 2 k gi 2 ^j + 24kL 2 gi Xzr Q 5 (z - ^ J - 36k 2 L 4 g lZ fir 3 (z - ^ + 24zfiL 6 k 3 

(r (Loh^g^o 8 - 9zr 7 gi - zr Q 6 (L z 2 + 2L a z + 3/i + 6A - 9) 
-l8z\L 2 k gi r 5 + 6zL 2 kr 4 + 27 L 4 k 2 fizr 3 gi + ftzfiL 6 ^))" 1 
.92 = (2. 9 ir 12 (L (gi + h x 2 (z 2 + l) L ) 

+z (-3 9l 2 X - 2h 1 2 g 1 L - 2h 1 4 L 2 )) - 3z 5l r n (z (2 9l + SL^ 2 ) - 13 5l + hh x 2 L ) + 

+ r 10 (-2z 4 L ( gi + h^Lo) + 2z 3 L (2 9l - h^Lo) 

+ z 2 (ALo^kXh^gi 2 + gi (42/i + 20A + 2) + 2h 1 2 L a (3/i - 2A + 7)) 

+z (6 gi 3 L 2 k (3/i - 2A 2 ) - StfkXh^g^Lo + 9l (24/< + 48A - 72) + 6/ii 2 i (/* + 2A - 3)) + 4L 2 fcA/ii 2 . 9 i 2 i ) 

- 6L 2 fcz 5 iAr 9 (z (Agi + bh 2 L a ) - 22 5i + bh 2 ) 

- 4L 2 kr 8 (z 4 XgiL - 2z 3 X gi L + z 2 {^iL 2 kh^ 9i 2 Lq + g x (-3 - A - 10A 2 - 21A/i) - 3/ii 2 L ^ 

27 '-2„3\,, or_„r2i. 2„ 2 , o„ / i2 t,,\ , 10 \ , q\ , o l 2 t- \ . 3 r _ , , j 2 . , . 2 „ 2 



+z ( -—kL 2 gi 3 X^ - 3L [iL 2 khi 9l 2 + 3.9i (-8A 2 - 4/tA + 12A + 3) + 3h 2 L j + -L^L 2 kh 2 gx 



+ 9L 4 zk 2 gi r 7 i^z i gi Ufi - -A 2 J + 5/i/ii 2 ioJ + 9i (-18/i + 12A 2 ) - 5/i/ii 2 £o 

+ 6L 4 zk 2 gi r Q e (z 3 fiL ~ 2z 2 ^L + z (-21/i 2 -fi- lOXfi + 4A) - 9kL 2 ^ gi 2 - 12/i 2 + /i (36 - 24A) - 12A) 
+ 36/iL 6 zfc 3 3l 2 Ar 5 (2z - 7) - 12^iL 6 zk 3 r Q 4 (z {2g x - L /ii 2 ) - 8.91 + /ii 2 L ) 
+27 fi 2 L s zk 4 gi 2 r 3 (2z - 5) + 2Az gi ^L 8 k 4 Xr 2 (z - 1) - ZQzg^ 2 L 1 ^ {z - 1)) 
• (r (-r 4 - 2AL 2 fcr 2 + 3L 4 k 2 ^) (-L o9l (z - 1) h^r^ + 9zgW 

+ zr a 6 (z 2 L Q + 2zL Q + 3/i - 9 + 6A) + 18zXL 2 k gi r 5 - 6zL 2 kr a 4 - 21L 4 k 2 ^zr 3 gi - Qz^k 3 ))^ 1 
h 2 = -h (r 5 (z {-2g x - h^Lo) - 3 9l + h^Lo) + 3zr 4 - 4L 2 k 9l Xr 3 (2z + 3) 
+6zXL 2 kr Q 2 + 3fiL 4 k 2 gi r {2z + 3) - 9L 4 k 2 [iz) 
■ (r 2 gi (-2r 4 - 4AL 2 fcr 2 + 6L 4 fc 2 M )) _1 



where g\ is given by (|23p . Each coefficient depends on the independent parameters r$ (the horizon radius) and h\ 
(proportional to the field strength at the horizon). 

Appendix C: Large r Series Solutions 

For large distances away from the black hole, r » L, we present here series solutions for 5 dimensions. 
Considering first the ansatz in (IM|) for k = 0, the field equations to first order of e imply 

= 2r 2 h e " + 2rhJ (z + 4) + zr (gj - fj) + 6zg e 

= 2r {z - 1) h e ' + 3rgJ + (z 2 - z + 12) g e + (z + 3) {z - 1) (2h e - f e ) 

= 2r (z - 1) h e ' + 3r.fJ + (z 2 + ^±^-l z - 18/x - 6^ g e + (z - 3) (z - 1) (2h e - f e ) 
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The solutions for f e {r),g e {r), and h e (r) yield integer powers of r in a number of special cases. For h e (r), we find 

(C 2 r~ 7/2 + C 3 r^ 2 ) (CI) 

/a (p 2 r-^ 2 + 2? 3 r 7/2 ) (C2) 

/ 2 (j- 2 r-^ 2 + J- 3 r-^ 2 ) (C3) 

where C%, C2, and C3 are integration constants, Lo is defined as in Appendix iBl and 

7 2 = { (-21/1 + 2A - 9) z 2 + (18/t - 20A + 26) z + 51/t + 50A - 33} (-1 + 3/t + 2A) _1 
-n ^ 2L (z - 3) (z - 1) 



MO 


= dr -3 -* - 


hr (-3- 




= Oir- 3 " z - 






= v ir - 3 - z - 





2L z 2 - (1 + 6A + 21/t) z - (9 + 6A + 45/t) 
P 2 = C 2 (£ 1 +£ 2 ) 
2? 3 = c 3 (A - £ 2 ) 

^2 = C 2 (£ 3 + £4) 

•7"3 = C3 (£3 — £4) 

£1 = ((1 - 8A - 21/t) z + 2 + 2A + 12/t) (• (36A 2 + (132// - 28) A + 9 - 30/t + 153/i 2 ) z 2 

+ (-40A 2 + (-72/t + 56) A + 12/i - 90/i 2 - 26) z + 68A 2 + (-92 + 132/t) A - 78/t + 33 + 81/t 2 ) 1/2 
£2 = (-8A 2 + (-10 - 78/t) A - 153/i 2 - 24/t + l) z 2 + 

(20A 2 + (132/t + 4) A + 5 + 261/i 2 + 42/t) z + 6 (6A + 16/t - 1) (A - 1) 
£ 3 = (z - 1) • £ • (■ (36A 2 + (132/t - 28) A + 9 - 30/t + 153/i 2 ) z 2 

+ (-40A 2 + (-72/t + 56) A + 12/t - 90/i 2 - 26) z + 68A 2 + (-92 + 132/t) A -78/t + 33 + 81/t 2 ) 1/2 
£4 = (z - 1) • L ■ ((-1 - 18A - 21/t) - 9/t - 10A + 27) 
1 



K 



2z (12z/t - z + 5zA + A - 2 - 3/t) £ 



For k = ±1, we can represent the asymptote functions as series 

2(n+z)-3 

/w = i+ E 7 

i=l 
2(n+ 2 )-3 

SW = 1+ E ^ 

2(n+z)-3 

Mr) = i+ E 5- 

1=1 

The coefficients can be determined from direct calculation. Due to the equivalence of our field equations, they match 
the values obtained in third order Lovelock gravity (l8| . once the substitutions /t£ 4 = — Q3 and A£ 2 = a 2 are made. 

For z = 2, we obtain nonzero coefficients only for the powers r~ 2 , r~ 5 , r~ 7 , r~ 9 . Just as in Lovelock gravity, at 
z = 2, all of the even powers of r until r~ 5 are present in the large r expansion of the asymptote functions. 
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